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Abstract 

The theory of neutrino flavor rotations is discussed in terms of wave function 

solutions to the Dirac equation with a neutrino mass matrix. We give a critical 

review of the nature of neutrino oscillations. 
PACS numbers: 1234.f 
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I. INTRODUCTION 



There is considerable interest in quantum interference patterns in space and time which 
result from the notion of flavor rotations. These include Kaon oscillations and B-meson 
oscillations which are thought to arise from quark flavor rotations. Some believe that Neu- 
trino oscillations also exist [1-4], but presently these are merely theoretical. Such neutrino 
oscillations are thought to arise from a possible flavor rotated Neutrino mass matrix. In all 
such cases, the quantum interference patterns are in part due to mass splitting, but they 
are also in part due to diffraction effects which would be present even if the mass splitting 
were not present. 

One should not forget that the invention of quantum mechanics was needed to explain 
amplitude interference oscillations without any recourse to the notion of mass splitting. For 
example, free electrons are described by a Dirac wave equation 

{-ih-f^d^ + mc) ipix) = 0, (1) 

where m is the electron mass. Using the usual spinors 

(7% + T^c) u{p) = 0, u\p)u{p) = 1, = -(?7^c)^ (1) 

one may construct a solution which consists of two plane waves 

V'(x) = Aiu{pi)e'P'-^/'' + A2^^(p2)e'^^•"/^ (3) 

and which yields an oscillating current 

j''{x) = ec4i{x)^'^'il){x) . (4) 

It is 

f{x) = ec (|yli|2^l(pi)7'^«(pi) + |A2|Mp2)7^«(P2)) + 2ecKe (yi;yli«(p2)7^«(Pi)e^''") , (5) 

where hk = pi~p2- Note that the space-time phase interference factor described by exp{ik-x) 
is Lorentz invariant, and that electron interference patterns in space and time are simply 
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and covariantly described by the current four vector j^{x). Although low energy electron 
diffraction (LEED) experiments carried out every day[5,6] are most often described in the 
laboratory rest frame, the description of the observed space-time oscillation in the current 
j^{x) is not overly difficult to describe in any Lorentz frame. Finally, the description of 
electron interference by normalized wave packets (explicitly including the spin s), 

i^{x) = J2[ A{p, s)u{p, s)e'^-^'HT, dT = , ^"^'^ u^p, s')u{p, s) = 5,s, (6) 

s V |Pr + ^^c^ 

by no means eliminates the oscillations in the current 

j'^(x) = ecJ2j J A*{p', s')A{p, s)u{p', s)-^^u{p, s)e'^P-P'^-^'^drdr. (7) 

s's 

While a sufficiently wide energy distribution in the wave packet will produce a low oscillation 
visibility, LEED machines routinely produce an electron beam which allows for the easy 
observation of electron diffraction via the factor exp(i(p — p') ■ x/U) in Eq.(7). 

We have reviewed these well known features of electron diffraction oscillations mainly 
because the case of massive neutrinos has for the most part (in the literature) been treated 
by a completely different set of rules than those implied by the Dirac equation. Since we 
feel strongly that the Dirac equation is perfectly adequate to the task of describing freely 
moving spin 1/2 particles, including massive neutrinos, we wish to compare the Dirac theory 
to those other theories which rely on a more obscure formalism. 



II. DIRAC EQUATION 



In the Dirac theory, the massive neutrino wave functions may be denoted by na{x) where 
a — e, ii,T denotes the flavor index and nia are the neutrino masses. The Dirac equation 
reads 



n^{x) 
nT-ix) 



\ ( 

+ 



mgC 



m^c 



ne{x) 



\ 



m^-c 



ni^[x) 



nT-{x) 



0. 



(8) 
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where the wave function has twelve components, i.e. four spinor components times three 
flavor components. The physical neutrino wave function ui^x) of three possible flavors is 
produced by the rotation 



^ee ^efj, -^er 



ne{x) 



R^e R^fi RfiT 
Rq-e RtiJj Rtt 



and obeys 



{-%h-i^d^ + Mc) v{x) = 0, 
where the neutrino mass matrix obeys 

Mah= RacrricRcb: o, 6 = (e, t) , 



c=e,n,T 



and the matrix R is unitary 



(9) 



(10) 



(11) 



(12) 



The general solution of the massive neutrino Dirac equation is best discussed in terms of 
the propagator. 

III. NEUTRINO PROPAGATOR 



The neutrino propagator is a twelve by twelve matrix obeying 

-tl^d^ + ^]G{x- x') = 5{x - x'), 



h 



or equally well 



(id'^Gix - a;')7'^ + G{x - ^0^) = S{x - x'). 
Consider a space-time region Q. Prom the obvious identity 

u(x) — / 5(x — x')iy(x')d'^x' , x E 
Jn 

and Eq.(14) one finds 

iy(x) = (id'^G{x - x')Y + G{x - x')^^ u{x')d^x\ xeVL, 



(13) 



(14) 



(15) 



(16) 



and 

u{x) = [{id'^G{x - x')Y) i^ix') + G{x - x') {iYd'^iy{x'))] rfV, xefl, (17) 
where Eq.(lO) has been employed. Eq.(17) imphes 

u{x) J^d'^{G{x -x')-f^u{x'))d'^x', x e n, (18) 

which may be converted into a "three surface" integral on the boundary dfl of the region 
fl, i.e. 

u(x) = i £^ G{x - x')-f^u(x')d^T,'^, x e Q. (19) 

From Eq.(19) it is evident that the propagator allows one to compute the full neutrino wave 
function for all x e if the wave function is known on the boundary x' e Writing 
Eq.(19) with flavor indices made explicit yields 

iya{x)=i Yl { Gab{x - x')j''ub{x')d%, X G Q, a = (e, t) , (20) 

b=e,iJ,,T 

where 

Gab{x -x')= Yl RacRcbS{x - x' ] TTlc), (21) 
c=e,iJ,,T 

and S{x — x';m) is the ordinary Dirac-Feynman propagator 

(^-i^^d^ + S{x -x';m)^S{x-x'). (22) 
The complete solution to the neutrino wave function problem is then formally 

^a{x)=t EE/ RacS{x-x'-mc)YR;b''b{x')d%, X e Q. (23) 
Eq.(13) can be solved[7] in the form 

G{x-x')^ (i^d, + ^^V{x-x'), (24) 

where 

(-9^9^ + {Mc/hf) V{x - x') = 5{x - x'). (25) 
From Eq.(19) and (24) it follows that 

u{x)=(trd, + ^^ip{x), (26) 



where 

(p{x) = i 



£^V{x - x')Y'^{x')d^^'i^, xen. (27) 



Eq.(27) is a proper starting point for discussing possible neutrino oscillation experiments. 

IV. OUTGOING NEUTRINO WAVES 

Here we consider neutrinos (going forward in time) in a fixed Lorentz frame. Anti- 
neutrinos (going backward in time) may be discussed similarly. Let i/(r, 0) be the initial 
neutrino wave-packet function at time zero in the Lorentz frame of interest. At later times, 
Eq.(27) reads (with (3 = 7°) 

(p{r, t) = i jvir- r', t)/3i/(r', 0)dV, (t > 0). (28) 

From Eq.(25) 

^ ' ^ i-oo V 47r|r-r'| J \2nhcJ ^ ' 

where the flavor matrix 

cP{E) = sJe-' - (Mc2)2. (30) 



Putting 



one finds 



^ /_>.(r)e-/» (^) , (31) 



' ApiV(E)\r-r'\/n' 



M^) = I [ ^ . /5Kr', 0)d'r'. (32) 



For large distances 



Mr)-[—^jf{E,f), f=^, |r|»|r'|, (33) 
with the production amplitude (twelve components for spin and flavor) 

/(^,f) = J e-'^(^)^-^'/?i^j,(r',o)(i3r'. (34) 
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In space and time, the asymptotic outgoing wave is 
Employing Eq.(ll), one finds with flavor indices made explicit 

I I o=e,/i,T c— e,//,T \ / 

where 

cp,{E) = yjE^ - {m,c^y. (37) 

Evaluating the energy integrals on the right hand side of Eq.(36) by steepest descents, one 
finds the energy 



dE 

yielding 



^ME)|r|-Et)^_ =0, (38) 



^ = Vb, Ef, = , Pb{Eb) = J 39) 

t ^1 - {vb/cY ywV^ 



so with the proper time defined by 



'piXEi)\A - Ebt' 



Tb = t^l - [Vb/cf = - ^-^^^^ ^ , (40) 

^ \ &mb ) 

Eq.(37) reads (putting the integration variable E — Ei, + e), 



M^,t)^^ Y: E Rabe-''''"''^^/''Rbc'Fc(Eb,t-^-^,r] 

1^1 b=e,n,T c=e,iJ.,T \ Vb J 



(41) 



where 



F.W,a) = £.-«/»m + ..f)(5^). (42) 

Eq.(41) is the central result of this work. 

Eq.(41) contains the phase factors ex.p{—ic'^mbTb/h), which are determined by the phase 
velocities, and the group envelope Fc{Eb,t — (|r|/i>6), r)which propagates with the group 
velocity. This decomposition (phase velocity in the phase and group velocity in the enve- 
lope) has been standard quantum mechanics for well over half a century. See for example 
the standard scattering theory treatise of Goldberger and Watson[8]. We stress this point 
because in the literature some workers put the group velocity into the phase, thus forgetting 
why the phrase "phase velocity" exists. 



V. EXPERIMENTAL IMPLICATIONS 



In experiments, there is a total outgoing amplitude (for a neutrino of flavor a, when the 
initial neutrino had flavor c) given by Eq.(41); It is 



total 



(43) 



The absolute value squared of the amphtude is then 



otal 
c—*a 



b=e,^,T d=e,iJ,,T 

with the neutrino oscillation phase factors 



Vb 



r\FjEd,t-—,r 



\r\ 

Vd 



(44) 



(pbd = (mdC^Td - mb(^Tb)/h. 



(45) 



The neutrino kinematics in the (laboratory frame) regime rribC^ ^ Eb must now be 
discussed. The velocities appearing in Eq.(43) are given by 



Vb 



A 



1 - 



mld^ 



c 1 



2 4 

2£f 



+ ... , {mbC^<^Eb), 



(46) 



very slightly less than light velocity. For example, if (mbC^/e) < Ivolt, and (Eb/e) > lO^volt, 
then [(c — Vb)/c\ < 5 x 10"^^; i.e. the velocity (laboratory frame) is less than light velocity 
only by the thirteenth significant figure. In the CHORUS and NOMAD experiments where 
neutrino energies are much higher (~ Gev) than the above estimates, it is safe to ignore the 
differences in the wave-packets Fc in Eq.(44) and use the phase coherent flavor conversion 
transition probability 



'Pcoherenti.C ^ o) — E E ^abRlb^ld^cd 
b=e,ix,T d=e,ix,T 



(47) 



where Eq.(45) has been employed. Note that the use of many proper times in Eq.(47) (i.e. 
the fact that each of the neutrinos with mass mb,d has its own internal proper time Tb^^), 
makes the phase coherent flavor conversion probability in Eq.(47) Lorentz invariant. 

This coherent flavor conversion probability is not so obvious when neutrinos come from 
distant stars. Intuitively, one would expect that the neutrino wave packet of a heavy massive 
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neutrino arrive on earth later than the wave packet from a hght neutrino, if the original 
neutrino source was (say) an exploding star. The wave packet peaks for two different mass 
neutrinos would be separated in space by Ar ~ r(Av/c) where Av is the difference in the 
two neutrino velocities and r is the distance from the distant star to the earth. This means 
that 



Ar ~ cm 



light — years 



Av 



X 10^^ {stellar .source). 



(48) 



Thus the phase coherence in Eq. (47) would be scrambled and one would expect (when the 
source is a distant star) an incoherent flavor conversion transition probability 



Vincoherent{c ^ a) = ^ \Rab\'^\Rcb\'^ , (stellar.source). 

b=e,n,T 



(49) 



Here Ar is so very much larger than the neutrino detecting nuclear event length scale in the 
target on earth. 

If the source of neutrinos is the sun, then the wave packet peak separation is given by 

-Av^ 



Ar ~ cm 



X 10^^ 



(50) 



■^earth— sun _ 

where Rearth-sun IS the distance from the earth to the sun. Again the incoherent flavor 
conversion transition probability applies; 



Vincoherentic ^ o) = ^ \Rab\^\Rcb\'^ , {sUU SOUrCc) . 
b=e,iJ,,T 



(45) 



Let us now compare the central result of this work, i.e. Eq.(47), to similar results as they 
appear in the literature[9-ll]. 

VI. PREVIOUS ABUSES OF DIRAC NOTATION 

The conventional textbook discussions of neutrino oscillations with flavor mixing are 
simply described as follows: In Dirac notation, the three flavors of neutrino {u^, i^^, Vt) are 
related to those neutrinos of fixed mass 



M\na >^ Ma\na > , {a^e,ii,T), 



(46) 
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by a flavor rotation matrix via 

|i^a >= Rab\nb>, \na>= iR~^)abWb>, {a^e,ii,T). (47) 

b=e,fj,,T b=e,ix,T 

In the course of time, the state of a given flavor neutrino changes via 

exp (-iiJi/^) |i/c >= Y RcbG:xp{-^Ebt/h)\nb >, (48) 

b=e,fj,,T 

where 

= V^M + c'IpI', (fe = e,//,r). (49) 
and where p is the neutrino three- momentum. Eqs. (47) and (48) imply that 

< i/„|exp(-im/^) >= Y Rcb{R~\aew{-iEb{p)t/h). (50) 

b=e,fi,T 

One then concludes that the transition probabilities between neutrinos of different flavors, 

Pwrong{c -^a,t) = \<ua\ exp {-iHt/h) \u, > f ????, (51) 

oscillate in time. The above described conventional arguments are clear, elegant, appealing, 
and wrong. 

One may suspect that the conventional discussion above has some problems with Lorentz 
invariance by applying it (as very many workers have done) to a neutrino produced in the 
sun and detected on earth. Suppose that a nucleus "at rest" in the sun flres off a beta-decay 
neutrino. There will be according to the above (incorrect but conventional) discussion a 
neutrino superposition of three different energies Ef,, {b = e, /i, r) but all with the same 
three-momenta p. Not all of the nuclei on the sun are "at rest" . Now suppose that a moving 
nucleus flres off a beta-decay neutrino. Can one believe that the three energies are yet again 
all different but the three-momenta are all yet again the same? The answer is obviously no! 
If Lorentz symmetry is invoked, then momentum and energy together form a four vector 
p — {p,E/c). One should put — —{rribcY for an on mass shell neutrino, and which of 
the four momentum components may be the same or may be different will then depend 
on the reference frame. At least the possibility that the neutrino is in a superposition of 
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different three-momenta should be considered. This requires that Eq.(49) for an on mass 
shell neutrino should be replaced by 



With this replacement the conventional discussion is raou; clear, elegant, appealing, and still 
wrong. 

Exactly what is the amplitude < Vfinail ^^p{—iHt /h)\i'initiai > supposed to mean? If we 
assume that the neutrino is a simple spin 1/2 particle without further internal structure 
(beyond the flavor label), then perhaps it is supposed to mean something like 



< ly final] exp{-iHt /h)\Uinitial >= i y]inal{r) ^W{-iHt /fl)PinUial{r)d^r ???? (53) 



We have chosen a time t (which means we have chosen a Lorentz reference frame), and a 
standard Dirac Hamitonian (with flavor indices implicit for example in the mass matrix Ai) 



If the above is indeed the case, then there is a remaining problem. In Eq.(53), the inner 
product contains an integral over space. 

Space is not the same as time. When there are energy superpositions there will be 
time oscillations. When there are momentum superpositions there will be space oscillations. 
Furthermore, when the neutrino comes from the sun to the earth, we know mainly where 
the neutrino was detected. The predicted oscillations should be (at least in part) in space. 
An integration over space seems a peculiar way to get quantum interference (oscillation) 
diffraction patterns that exist in space. Consider electron diffraction (in real experiments 
scattering electrons off a crystal face). If "thought experiments" are more appealing, then 
consider two-slit electron diffraction. To get the interference (oscillation) pattern in space, 
one absolute squares a spatial amplitude. If you integrate over space, then you get a to- 
tal probability of unity (which is true but not really of much use in discussing diffraction 
oscillations) . 




(52) 



(54) 
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Our suggestion has been in the previous sections the following: (i) If one wants to study 
neutrino oscillations in space-time, then it is best to study a wave function which actually 
depends on space and time; e.g. try the Dirac equation Eq.(lO). (ii) If one wants to study 
how the neutrinos are distributed in space and time, then it is best to study the current 
j'^{x) — cU{x)'y'^h'{x) which actually describes the distribution in space and time, (iii) Finally, 
if one wants to study how a neutrino can propagate from the sun to the earth (or across a 
large laboratory), then it is best to study the propagator Eqs.(13) and (19) which actually 
describes how the neutrino propagates. 



VII. CONCLUSION 

The coherent neutrino oscillation flavor conversion probability 

n J ^^^^ 

depends on the proper times of the massive neutrinos; i.e. 

%(t)b = -mbC^Tb = {pb\r\- Ebt). (56) 

Some previous workers employ only the energy-time factor Ebt in the phase 4>b, while other 
previous workers employ only the momentum-space factor pb\r\ in the phase 4>b- To consider 
only energy-time or momentum-space and (not both) tells only half the story. The true 
Lorentz invariant phase 0;, requires both. 

Finally, neutrino flavor conversions may be described sometimes by Vcoherent{c a) 
and sometimes by Vincoherent{c — a), depending upon what is assumed about neutrino 
masses, rotation matrix elements, energies of neutrino sources and length scales of actual 
experiments. There is no universal simple formula which holds true in all regimes. 
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